Small-angle impurity scattering and the spin Hall conductivity in 2D systems 
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An arbitrarily small concentration of impurities can affect the spin Hall conductivity in a two- 
dimensional semiconductor system. We develop a Boltzmann-like equation that can be used for 
impurity scattering with arbitrary angular dependence, and for arbitrary angular dependence of the 
spin-orbit field b(k) around the Fermi surface. For a model applicable to a 2D hole system in GaAs, 
if the impurity scattering is not isotropic, we find that the spin Hall conductivity depends on the 
derivative of b with respect to the energy and on deviations from a parabolic band structure, as well 
as on the angular dependence of the scattering. In principle, the resulting spin Hall conductivity 
can be larger or smaller than the "intrinsic value", and can have opposite sign. In the limit of small 
angle scattering, in a model appropriate for small hole concentrations, where the band is parabolic 
and b oc k^, the spin Hall conductivity has opposite sign from the intrinsic value, and has larger 
magnitude. Our analysis assumes that the spin-orbit splitting b and the transport scattering rate 
r^^ are both small compared to the Fermi energy, but the method is valid for for arbitrary value of 
6r. 

PACS numbers: 72.25.-b, 73.23.-b, 73.50.Bk 
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I. INTRODUCTION 

An electric current passing through a semiconductor 
can induce spin polarization near lateral edges of the 
sample, with opposite sign at opposite edges. Such an 
effect, commonly referred to as the spin Hall effect, had 
been predicted more than 30 years agc^ and was ob- 
served recentlyiSi2i4 This effect results from the coupling 
between spin and momentum of an itinerant electron. 
Two specific mechanisms have been considered for the 
spin-Hall effect. The "extrinsic" mechanism couples spin 
to momentum during events of impurity scattering in 
the Mott skew scattering channeli*^. As a result, scat- 
tered electrons of different spin polarizations propagate 
towards opposite boundaries of the system. As opposed 
to the extrinsic, the "intrinsic" effect originates from a 
spin-split band structurei^ rather than from impurity 
scattering. 

Experiments have been performed on both electron- 
dopedi2i4 and hole-doped^ GaAs-based semiconductor 
structures of different dimensionalities. A spin Hall ef- 
fect in three-dimensional n-doped GaAs films^ has been 
explained in terms of extrinsic skew and side-jump impu- 
rity scattering^'"'^^ while the intrinsic mechanism was es- 
timatedii to yield spin accumulation smaller by an order 
of magnitude than actually observed in the experiment. 
The two-dimensional electron gas on a (110) surface of 
GaAs, studied in Q, was also in the dirty hmit, and the 
spin Hall effect observed there was attributed to extrin- 
sic effects. On the other hand, p-doped two-dimensional 
(2D) samples used in the experiments of Ref. Q were 
estimated to contain very few impurities and it was sug- 
gested that the observations reflect a spin-Hall effect of 
the intrinsic type. 

While experiments directly measure spin accumulation 



along the sample edges, theoretical contributions to the 
field of spin-Hall effect are mostly concerned with a dif- 
ferent quantity which is easier to calculate, namely, the 
spin current. The latter is conventionally defined as the 
expectation value of the operator — C{si,Vk}, where 
s and V are, respectively, the operators of the carrier spin 
and velocity (which is, in general, also spin-dependent), 
and the normalization constant C has been chosen differ- 
ently by different authors. (In the present paper, we de- 
fine Si to have eigenvalues ±?i/2, and we choose C = 1/2, 
so that the spin current is equal to U/2 times the differ- 
ence in the particle currents for carriers with opposite 
spins.) Apart from the normalization, the definition of 
spin current is largely a matter of convenience since it 
is not rigorously related to spin polarization (accumu- 
lation) by a conservation law, as opposed to e.g. the 
electric charge conservation, namely, dtSi + ^kjl 0- 
Non-conservation of spin-current (due to spin precession) 
makes the problem of its conversion into polarization 
highly non-trivial. Moreover, since spin current is even 
with respect to time inversion, it is not necessarily absent 
even in equilibrium-—. In the present paper we concen- 
trate on the 2D electron and hole systems. In this case 
the components of spin current polarized along the 
direction (z) perpendicular to the 2D plane appear due 
to non-equilibrium conditions only. 

The paper ,8] predicted a universal value for the intrin- 
sic spin-Hall current in a 2D electron gas, with spin-orbit 
coupling of the Rashba form, which, using the normal- 
ization of the present paper is given by 
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independent of the strength of spin-orbit interaction. 
(The prediction of Ref. [2| for three-dimensional hole 
systems, though dependent on the electron density, was 
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also independent of the strength of SO couphng). These 
results immediately posed a question of how impurities 
would affect the spin-Hall conductivity. The derivation 
of formula neglected scattering, but also assumed a 
steady state. The latter, however, is impossible in an 
infinite system if momentum relaxation is absent. Soon 
a number of papers appeared about the role of impurity 
scattering, which, initially, reached differing conclusions. 
Some studies encountered decreasing of the magnitude 
of the effect with increasing scatteringiii4, while others 
found a universal value ^ even in the presence of scat- 
tering. However, it is now generally agreed that in the 
simplest model of a 2DEG, with a spin-orbit coupling 
that is linear in the wavevector k, the dc spin Hall con- 
ductivity will vanish, even in the case of arbitrarily small 
impurity concentration. Inoue et al^^ first reached this 
conclusion in the clean limit, 6t — > oo, where the spin or- 
bit splitting b exceeds the transport scattering rate t^^, 
while entertaining the possibility of a finite effect for fi- 
nite values of br. 

The present authors!^ developed a theory based on 
a quantum kinetic equation to describe spin-polarized 
transport in the two-dimensional electron systems for ar- 
bitrary ratio of spin-orbit splitting and scattering rate, as 
long as they are both small compared to the Fermi energy 
{b,T~^ Ep). Applying this theory to the spin-Hall 
effect in a 2DEG with (k-linear) Rashba coupling and 
isotropic impurity scattering, we established the vanish- 
ing of spin Hall conductivity for arbitrary value of the 
product br. However, the cancellation is complete only 
in the bulk. Near the contacts which inject unpolarizcd 
electrons, spin-currents are non-zero and, hence, a spin- 
polarization normal to plane can accumulate near the 
corners of a sample. 

The models used in Refs. ^| and [31 assumed short- 
range impurity potentials, with isotropic scattering. The 
result of vanishing spin Hall conductivity in a 2DEG with 
Rashba coupling was subsequently confirmed by several 
studiesiiiiSiiSiS^, and were shown to hold for arbitrary 
angular dependence of the impurity scattering. Early 
numerical calculationsSSiSi seemed to contradict the the- 
oretical work but later ones supported the vanishing of 
spin-Hall currcnt^^. Results for a 2DEG with Rashba 
coupling can also be applied directly to a model for a 
2DEG on an (001) surface with pure k-linear Dresselhaus 
coupling, as the two problems are related by a rotation 
in spin space. 

As was pointed out by Dimitrova^^, the vanishing of 
dc spin currents in an infinite 2D system with k-linear 
spin-orbit coupling has a simple explanation in terms 
of the following operator argument. Consider a system 
with the spin- independent disorder and the "Rashba" - 
type spin-orbit coupling^^ Hso — Oi{axky ~ (Tykx), where 
(T represents the set of three Pauli spin matrices, and 
k is the canonical momentum. The equation of motion 
for the operator of spin polarization yields {I — x.y): 
dtO'i = z[7?so, c;]/?! = —Aamji/h^, where m is the elec- 
tron effective mass. It follows from this identity that a 



nonzero expectation value of the spin current would re- 
sult in a time-dependent spin polarization along the cor- 
responding direction Z, which is impossible in a steady 
state. This operator argument applies also to a model 
with a combination of Rashba and k-linear Dresselhaus 
coupling. In agreement with this argument, analytical 
calculations confirm the vanishing of weak-localization 
(of the first order in l/Epr) corrections to the spin Hall 
conductivit}i2^. 

It should be emphasized that this vanishing of the spin 
Hall conductivity is characteristic of the 2D electron sys- 
tem with SO coupling linear in momentum. Cubic Dres- 
selhaus terms in the spin-orbit Hamiltonian, which are 
oc axkxky — aykyk^. for a 2D electron system on a (001) 
surface, should result in a finite spin Hall conductivity^-. 

For a 2D hole system in an asymmetric confining well, 
the Rashba effect gives rise to a spin orbit field b(k) 
whose direction winds three times around a circle in 
the x-y plane as k moves once around the Fermi circle. 
It turns out, in this case, that for short-range impuri- 
ties there is no cancellation in the clean limit, and one 
then obtains the full value of the intrinsic spin Hall ef- 
fect, at least in a simplified model which neglects cubic 
anisotropy.^^ Curiously, in contrast to the 2DEG, it is 
the vertex corrections that vanish in this case, analogous 
to the absence of the vertex corrections to the electrical 
conductivity for short-range impurities^. 

Vertex corrections are also absent in an n-doped 3D 
system with Dresselhaus interaction, Hso = ^^^xkxik"^ — 
kf ) -|-cycl. permut., and short-range impurity scattering, 
in this case due to the cubic crystal symmetry. In the 
absence of vertex corrections, the impurities lead to a 
smooth relaxational suppression of the spin-Hall conduc- 
tivity due to the broadening of the spectral function of 
carriers. In principle, one should obtain the full intrinsic 
value of the spin-Hall conductivity in the clean limit, but 
3D systems are always in the dirty limit, where the intrin- 
sic spin Hall conductivity is reduced by a factor ~ {br}'^- 
As mentioned above, the spin Hall effect observed in 3D 
n-type GaAs is most probably due to extrinsic scattering 
effects, omitted from the model discussed herei2i2ii£ 

It is natural to ask whether the absence of vertex 
corrections for the 2D hole system is special to the as- 
sumption of short range impurities, or whether it applies 
equally well to a model with arbitrary angular depen- 
dence of the impurity scattering. In the present paper, 
we develop equations to calculate the spin Hall conduc- 
tivity for arbitrary form of the impurity potential, and 
for arbitrary form of the spin-orbit field b(k) in a 2D 
system. Solution of these equations is simple in models 
which have an overall circular symmetry, and we obtain 
analytic results in the limit of small angle scattering, re- 
sulting from smooth disorder. This type of disorder is 
a good assumption for many 2D systems in which the 
impurity centers (dopants) reside relatively far from the 
heterostructure interface. 

The result of our analysis is that the spin-Hall conduc- 
tivity, in general, does depend on the form of the impurity 
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scattering, and except for the case of isotropic scattering, 
depends also on the energy-dependence of the spin-orbit 
couphng and on the deviation from parabohcity of the 
energy dispersion for the carriers. In our circularly sym- 
metric model, the magnitude of the spin-orbit field b(k) 
is independent of the direction of k, while its direction 
winds N times around the unit circle in the x-y plane, at 
a uniform rate, as k moves around the Fermi circle. {N 
must be an odd integer.) The energy-dependence of b and 
deviations of the energy dispersion from a parabolic 
form near the Fermi energy, are then characterized by 
two parameters N and C, defined in Eq. (0J below, and 
for a given form of the impurity scattering, the spin Hall 
conductivity is found to depend on N, N and C- 

In Section II, below, we define precisely the models we 
are considering and the relevant parameters. In Section 
III, the transport quantities are expressed in terms of the 
Wigner distribution function, for which Boltzmann-like 
equations of motion are presented in section IV. In Sec- 
tion V, these equations are solved in the case of our circu- 
larly symmetric model. In Section VI, we obtain explicit 
formulas, Eqs. H57|l and (|51|l . for the two limiting cases 
of isotropic scattering and small angle scattering by the 
impurities. The relation of our results to previous theo- 
retical results and to recent experiments is discussed in 
Section VII. Some details of the derivations are presented 
in Appendix A, and an alternate derivation of the kinetic 
equations of motion, using Green's function methods, is 
presented in Appendix B. 



II. MODEL 

We consider a two-dimensional system of non- 
interacting electrons or holes, described by a one-body 
Hamiltonian of the form: 



H = ek 



1. 



;b(k) 



V{r) 



(2) 



where V is the potential due to impurities, and k is the 
kinetic momentum operator. In the presence of a vector 
potential A, we have k = p — eA/c, where p = — iV. 
Throughout this paper, e is the charge of the carriers, 
and we set h = 1 except where explicitly noted. We may 
include a uniform electric field E by letting A depend on 
time. 

By time reversal symmetry, the spin-orbit field b must 
satisfy b(k) = — b(— k) . Note that vectors such as k and 
r, which refer to spatial coordinates, have two cartesian 
components, while vectors in spin space, such as b have 
three components. 

We have assumed for simplicity that e^, the energy dis- 
persion in the absence of spin orbit coupling, is isotropic 
and a monotonic function of k. The formulas derived be- 
low can be extended in a straightforward manner to more 
complicated band structures, but the formulas would 
then be much less transparent, and the resulting integral 
equations are more difficult to solve. 



As a further simplification, we shall later specialize to 
models where b has a simple rotational symmetry; i.e., 
we assume 



0, bx + iby — bQ{k)e 



iNe 



(3) 



where 6o(fc) is a complex number, whose magnitude may 
vary with the magnitude of fc, but whose phase is inde- 
pendent of k. This will allow us to obtain analytic solu- 
tions to the equations. The final results will, of course, 
depend crucially on the winding number N . We shall 
see that results can also be affected by the energy depen- 
dence of b and by deviations of the energy dispersion 
from a parabolic form near the Fermi energy. We char- 
acterize these dependences by constants N and C, with 



N 



d\n\bo\ 
d\nk 



1 + C 



d\nv 
din A:' 



(4) 



for k equal to the Fermi momentum. In other words, we 
assume that near the Fermi energy, v oc and 



bx + iby oc fc^e' 



N9 



(5) 



For a 2D electron system on a (100) surface of a III- 
V semiconductor, the case = 1 corresponds to pure 
Rashba coupling, while A^ = — 1 corresponds to a k-linear 
Dresselhaus coupling. The case A^ = 3 arises in a circu- 
larly symmetric model of a 2D hole system, which ignores 
warping due to the tetragonal symmetry. In the limit of 
small hole doping, the band structure is parabolic and 
[bx + iby) oc {kx + ikyY] hence C = aiid N — N ^ 3. 
At larger values of the doping, when the Fermi-energy 
becomes comparable to the splitting between the light 
and heavy hole bands, the values of N is reduced, and 
the value of ( will be negative. At still higher doping, 
both light and heavy hole bands become occupied, and 
the situation is more complicated. 

The Hamiltonian Q assumes that there is no direct 
spin-orbit coupling associated with the impurity poten- 
tial or with applied electric field. For example, it omits 
terms of the form Act • (k x W), which are generally 
present in systems with spin orbit coupling. This term 
leads to skew-scattering by the impurities, giving an ex- 
trinsic contribution to the spin Hall conductivity. Since 
there is no direct general relation between the coupling 
constant A and the magnitude of the intrinsic spin-orbit 
splitting b, it is at least logically consistent to consider a 
model with A = 0. Moreover, the skew scattering cross- 
section is higher order in the impurity potential V than 
the ordinary transport scattering cross-section, so it is 
normally dominated by places where the carrier enters 
a region of strong potential gradients, close to the im- 
purity. In a remotely doped 2D hole system one might 
expect that the skew scattering contribution to the spin 
Hall conductivity should be relatively small, at least if 
one can ignore scattering by residual impurities close to 
the 2D system. 
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We concentrate on systems which are large compared 
to the mean free path. The "mesoscopic" spin-Hall 
eSec^^^t^^^ in ballistic heterostructures is, therefore, be- 
yond the scope of the present paper. 



III. BOLTZMANN FORMULATION 

We begin by defining a Wigner distribution function 
for the carriers, which is a 2 x 2 matrix in spin space, 
given by 

1 



s 



S T 



(6) 



In thermal equilibrium, if spin orbit coupling is absent 
(b — 0) and impurities are neglected, we have the stan- 
dard result 



the two spin states cr = ±1, which may be written as 
(27r)-2(fce + Sk)/{ve + Sv), where Sk = {ba/2ve) and 
Sv = {d'^ek/dk^)5k -f {l/2)adb/dk. The shifts in mo- 
mentum and velocity have the opposite sign for the two 
values of a, and the final term in Eq. (|10() is obtained by 
expanding to first order in these shifts. 

Finally, we introduce functions n{9,e) and ^{9,e), 
which describe the excess number and spin densities in 
momentum direction 9 and energy e, defined by 



naf}{9,t) - n''^p{9,e) 



n{9,e) 



6at3 + ^{9,e)-(Ta,p. (11) 



If n and # are known, we can readily compute the 
particle current and spin current, as well as the particle 
density and spin density at a spatial point r. We define 
the spin density and spin current by 



(12) 



a, 13 



nai3 = /o(e)'^(e - £k)Sal3, 



(7) 



where /g is the Fermi function. For 6^0, there are 
two different wave vectors at a given energy e and a 
given direction 9, corresponding to spin states parallel 
or antiparallel to b. Introducing an index cr = ±1 to 
distinguish these two states, one finds in equilibrium, 
«Q/3 = n^'^(fc,e) with 



(T = ±l 



(8) 

Impurity scattering broadens the S functions in these ex- 
pressions, but does not shift their centers or alter their 
amplitudes, provided that the scattering rate is small 
compared to the Fermi energy Ep. 

We define the distribution as a function of energy and 
direction in k-space by 



^a/3(^', e) = 



k dk 
(2^ 



Hapik, e). 



(9) 



with k = k{cos9, am9). (We suppress here the indices r 
and t, and we shall omit other indices as well when the 
meaning is clear.) Substituting (|HJ into we find, in 
equilibrium, to first order in h/Ep, 



F(r) = i^(V^^(r)K,v}„^^^(r)), (13) 



a, 13 

where v is the velocity operator, 

_ dek CT 9b(k) 



(14) 



ak 2 ak ■ 

At any point r, we then have for densities and currents 
(n) -no= d9 den{9,e), (15) 



/ d9 de$^(0,e). 



d9 



de v^{9)n{9,e), 



(16) 



(17) 



27T poo 1 

^ - / d9 de [M9)<f^{9, e) + -6v^{9)n{9, e)] 

J-oo ^ 



(18) 



Here 



(27r)2t;, 



■/o(e)<5a/3 



(10) 



where = and /cg is the value of k such that 
is equal to the given value of e. For e at the Fermi 
energy, fc^ is the Fermi momentum in the absence of 
spin-orbit coupling, and v,: is the corresponding Fermi 
velocity. The last term in (|10|l may be understood as 
arising from the difference in the densities of states for 



Ve(6') = We (cos 6*, sin 6*), 



is the "bare" velocity, and 
5v^i9) EE 



2 dk 



(19) 



(20) 



is the spin-orbit contribution. In the equation for j, we 
have dropped terms of order This is justified be- 

cause we assume |fo|/i?F ^ li Sind in the situations we will 
consider, values of |$| are smaller than the characteristic 
values of n{9, e) by a factor of |t'|/£'F or more. 
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IV. EQUATIONS OF MOTION 



We will be interested in the linear response of the sys- 
tem to an infinitesimal applied electric field E, which 
might, in general, depend on space as well as time. In 
developing the equations of motion, then, we need only 
consider terms which are of first order in $ and n, or 
first order in E; we may ignore terms of order En or E^. 

In the absence of spin-orbit coupling, the linearized 
equation of motion for n{6, e, r, t) may be written as 



dn 


dn 




dn 




dn 


'di " 


'di 


V 


" 'dt 


-\ 

E 


" 'dt 



where 



dn 
'di 



= -v^k ■ — n, 
ar 



is the standard advection term, 
2eEk cos 6 



dn 
'di 



-dh{e) 
de 



(21) 



(22) 



(23) 



describes acceleration of particles by the external field, 
and 



dn{9) 



dt 



27r 



de' K{9^9')[n{9')-n{e)], (24) 



is the collision integral. The scattering kernel K is given, 
in the Born approximation, by 



K{e-e') = Wiq)kJi2TTV,), 

q = 2/c,sin(|6'-6l'|/2), 
W{q) ^ {\V{ci)\'). 



(25) 



Here ^(q) is the Fourier component of the scattering 
potential at wavevector q, and we have assumed that 
its mean-square value is independent of the direction in 
space. In (j^SJ, we have assumed that the electric field 
E is oriented along the x-axis. The equations assume 
that spatial variations are slow on the scale of the Fermi 
wavelength, and they may need to be supplemented by 
boundary conditions derived from microscopic equations, 
in the case of a sharp interface or sample edge. 

The equations given above for dn/dt remain valid in 
the presence of spin-orbit coupling, to first order in b/E-p. 
The equations of motion for however, are affected in 
an essential way by a nonzero b. We find here 



(26) 



a* 


a* 


a* 




a* 




a* 




dt 


dt 




4 

b 


dt 


-\ 

E 


dt 


scat 



'dt 



-Vek ■ 

or 



(27) 



'di 



.-bx* + „|bx|)jL (28) 



'di 



Ee 



dfo 
de 



d9 



(b sin ( 



(29) 



a*(6') 



dt 



/ d9' K{e^e')[^{9')^^{9)] 

Jo 

d9' [m{9,9')n{e')-m{9',9)n{e)] . 

(30) 



Spin-orbit contribution to the scattering kernel can be 
divided into two different parts: 



M = M'^ + M'", 



where 



v,K{9-0') d 

4k, a^ 



kMke,0) 



(31) 



(32) 



originates from spin-orbit correction to the density of 
states, and 



AkeVe 



9'\ dK{e-9') 



d9 



(33) 

is due to the spin-dependence of the momentum transfer 
entering the scattering matrix element. In the case of 
small angle scattering, Eqs. and which are cor- 
rect to first order in b, require not only that b/Ep << 1, 
but also b < vpq, where q is the typical scattering mo- 
mentum. 

A more detailed derivation of these equations will be 
given in Appendices A and B. We give here just a brief in- 
dication of the origin of the various terms. We then shall 
explore some consequences of these equations in simple 
cases, which will shed additional light on their meaning 
and will provide some non-trivial consistency checks. 



A. Origins 

The contribution given by (|28|l may be understood as 
arising from precession of the spin, induced by the spin- 
orbit field b. The second term on the right hand side 
of this equation is zero in any case where the direction 
of b is determined by the direction of k, independent of 
the magnitude of fc, and we shall omit it in the following. 
The electric field term (|29|l refiects the action of E on the 
equilibrium distribution (jS)) , as explained in Appendix A. 

The scattering kernel M entering H3U|) gives a contri- 
bution to d^/dt proportional to n, which is present even 
when $ = 0. The term M'', given by t\'62\\ . arises because 
there are two different Fermi radii at a given energy e and 
a given direction 0, corresponding to spin states parallel 
or antiparallel to b, and the densities of final states are 
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generally different for particles scattered into these two 
spin states. The second term M™, given by (|33|) arises 
when the scattering matrix element depends on q, be- 
cause the momentum transfer actually depends not only 
on 9 — 9' ^ but also on the initial and final Fermi radii, 
and hence on the component of the spin in the directions 
of b(6') and h{e'). 



origin, but will have a spin direction that depends on r 
away from the origin. If one moves out along a constant 
direction f, the spin direction will precess about the di- 
rection b(k'), where k' is a vector of magnitude and 
direction f. The polarization $(r, 0, e) produced by such 
a scattering wave solution will give a time-independent 
solution of the equations motion derived above, at least 
within the Born approximation. 



B. Application to simple examples 

First consider a situation where E = and n{9, e) is 
independent of 9. The equations of motion for $ then 
have a time-independent solution where 



Akf de \ Vf 



(34) 



A change in electron density independent of 9 is just what 
one would find if one makes a small change in the chem- 
ical potential, by an amount Sfj, such that 



djo 
de 



(35) 



The change in the equilibrium distribution produced by 
a shift can be obtained from by replacing /o(e) by 
fo{e—S^). We see that this would change the polarization 
by an amount precisely equivalent to H34|l . as is required. 

As another example, we may consider a situation where 
there is no impurity scattering as well as no electric 
field. We may then construct spatially-varying time- 
independent solutions of the equations of motion of the 
form 



*(r, 61, e) = ab + c 
where b is 



g cos((5k ■ r) + {k X g) sin((5k ■ r) 



(36) 



wnere o is a unit vector parallel to h{k^^9)^ g is a unit 
vector perpendicular to 6, a and c are arbitrary constants, 
which may depend on 9 and e, and 



(5k = (cos 6*, sin ( 

2Vf 



(37) 



This steady state corresponds to a situation where we 
have populated eigenstates of the Hamiltonian with en- 
ergy e and wavevector-direction 9, whose spin at r = 
is polarized in direction ab + eg. The spin polarization 
precesses with a spatial-frequency Jk as one moves away 
from the origin. 

Finally, we may consider a situation where impurities 
exist only in a small region about the origin, whose radius 
is small compared to l/5k but large on the scale of the 
Fermi wavelength. We can then construct a scattering 
wave solution of the Hamiltonian, for which there is an 
incident plane wave of definite energy e and wavevector 
k, with polarization parallel or antiparallel to b(k). The 
scattered wave radiating from the origin will have spin 
parallel to the spin of the incident wave for r close to the 



V. CASE OF CIRCULAR SYMMETRY 

We now specialize to the case where b(k) has the sim- 
ple dependence on 9 given by l|2Jl. We wish to find the 
linear response to a uniform electric field, in an infinite 
homogeneous system. 

When the angular dependence of b is given by ((SJ, the 
equations can be solved by Fourier transform in 0, i.e., 
by expanding ^{9) in circular harmonics: 



m— — oo 



K{9-9') = 



(38) 



Since $2 is real, we have $f„ = (<l>l„j) . Also, since K 
is real and an even function of 9 — 9' , we see that Km is 
real and Km = K^m- Equations 1)28(1 and 129|l can then 
be written as 



dt 



2i 



dt 
d^L 



dt 
d<^m 



dt 



Eer]om . . dfo 

2 [""1,(^-1-1) ~ "m,(Af-l)J -Q^, (.-jyj 



Vo 



(40) 



The scattering contribution depends on n{9,e). For a 
dc electric field, one obtains the standard result, 



n{9,e) = Et 



ekf 



dfo 



cos ( 



27r2 V de 

where t is the transport scattering lifetime, given by 



(41) 



i.2 

T 



d9 K{9) sin^ 



(42) 
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We then find 



dt 



d^r. 



dt 



eE 



where 



7(±)=7o' + 7o"±A-, 



(43) 



-f) , (44) 



(45) 



Note that km are dimensionless constants that depend 
on the shape of the angular dependence of the impurity 
scattering, but not on its overall strength. 
Solution of these equations gives 



djo 
de 



2D 



(54) 



D = {t-'^/2) [|6o|V2(fcAr+i + fcAr_i) + 2kN+lkN^ 



(55) 

Note that the right-hand side of 154(1 is real, so that ^\ 
is pure imaginary. 

Integrating H54|l over the energy e, we find the spin 
currents j| = and jy = — CTsh-E', with a spin Hall con- 
ductivity given by 



7o^ = 2^77or(i^i-ifN)(iV-C), 



7o" = 2^'?or 



Kn+1 " Kn-1 



(46) 



, (47) 



(48) 



with 



K{9-6') = tan 



dK{e - 9') 



(49) 



fTsH = TTeVF 



2D 



(56) 



In the clean limit, |5|t 3> 1, one can replace the denom- 
inator D in Eq. by |6oP(A:jv+i + kN_i)/2. Since 
q:(±) (X \bo\, we see that CTsH is independent of r and inde- 
pendent of the magnitude of |6| in this limit. By contrast, 
in the dirty limit, when |6|t <C 1, we see that CTsH oc jforp 

Note that spin-Hall conductivity cancels at = 1 for 
an arbitrary scattering kernel K{9 — 9'). Indeed, accord- 
ing to Eq. 1461) . 7q is identically zero for = 1, while 
7^ and A"" given by Eqs. and (gTjl are equal. Thus, 
both 7''~-', and a(-) are zero, therefore the first term 
in the Eq. (|56|l vanish. Also, fco = by definition (see 
Eq. (53)), and the spin-Hall conductivity is identically 
zero. This result is in agreement with a general argu- 
ment given in Ref. p^. 



It is seen by inspection that the solution of these 
equations has all components zero except for $Ar±i and 
$1 = (<&li)*, which are therefore determined by a set 
of three coupled linear equations. For a dc applied field, 
where d^/dt = d^^ /dt = 0, these equations may be writ- 
ten in the form 



2i^l = {^b*o<pN+i + 6o*^-i)t, 



(50) 



Ec ( —0 f 



2 \ de 



Ee f-dfo\ , 



2"i^r(-) 



«(-), (52) 



where 



km = 2tt{Ko - Km)T, a(+) = [+7(+) +m{N+ 1) 



-7(-)+'7S(A^-l) -(53) 



VI. ISOTROPIC SCATTERING AND SMALL 
ANGLE SCATTERING 

For the case of isotropic impurity scattering, we have 
Kq — (27rr)~^, and Km = for m 0. We also have 
K ^0,j^^0, and 7(±) = 0. Then, for \N\ ^ 1, we find 
in the clean limit, 



CTsH 



eiV 
8^' 



(57) 



This is just the "universal intrinsic value" of the spin Hall 
conductivity for this model, as has been found previously 
by other authors, written for our present normalization. 



«(+), (51) The result will be reduced by a factor of (1 -I- |6| 



OI T 



if |6o|''' is not small. 

In the limit of small angle scattering, we have 



/•n" rp. 
K{9)-d9, 



K{9) 



m 

2^' 

9 d 
2d9 



K(9) 



(58) 
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FIG. 1: Spin-Hall conductivity CTsh, in units of e/Svr, for vari- 
ous values of the band curvature parameter i^, as a function of 
p, the fraction of transport scattering rate due to small-angle 
scatterers. All curves are for parameters N = N = 3. (See 
Eqs. (3) and (4) for definitions.) The clean limit, (bor ^ 1), 
is assumed. The case = corresponds to a circularly- 
symmetric model of a 2D hole gas in the limit of low carrier 
density. 



To compute Km-, one has to integrate by parts and ex- 
pand the integrand to the second order in 9. The result 
is: 



K„ 



We then find 



7o 



Ko 3m2 

— - H 

2 47rT 



%(iV2-l)(iV-C), 



(59) 



7o = 3?7o, 



(60) 



Using Eqs. gSJ, lO, (EHJl and wc thus obtain, in 
the limit of small angle scattering and \bo\T » 1, 



CTsH 



(61) 



As expected, this result gives CTsh — for |iV| = 1. 
For \N\ ^ 1, the result for small angle scattering can be 
larger or smaller than H57|) . and can even have opposite 
sign, depending on the values of N and C. 

For a 2D hole system in GaAs, at low doping, with N — 
N — 3 and C = 0, we obtain, for small angle scattering. 



3 e 

O-sH = -7-, 
5 TT 



(62) 



which has different sign and is approximately twice as 
large as the result l(57|l for isotropic scattering. 

In a remotely-doped 2D electron or hole system, the 
charged impurities that compensate the carriers are lo- 
cated far from the 2D system, so they contribute only 



long-wavelength Fourier components to the disorder po- 
tential, and produce only small angle scattering. As the 
set back distance is made large, their contribution to the 
transport scattering rate decreases, so that eventually 
the latter may be dominated by a small concentration 
of residual impurities close to the layer, whose scattering 
amplitudes may be nearly isotropic. Thus, it is natural 
to consider a model where both types of impurities are 
important, so the scattering kernel K{9) contains both 
a narrow peak at « 0, and an isotropic 6'-independent 
part. If the two contributions are weighted so that a frac- 
tion p of the transport scattering rate arises from small 
angle scatterers and a fraction 1 — p arises from large 
angle scatterers, then doing the Fourier transform, one 
finds that the coefficients km entering H56|) are given by 
km = 1 - p + pm?, for m 7^ 0. [From and (jHSJ, it 
may be seen that for arbitrary K{9), the transport scat- 
tering rate is determined by Kq — Ki , so that one always 
has fci = 1.] The constants 7oj7o'j and A™, given by 
Eqs. H46|) - (|48|l have no contribution from the isotropic 
part, and are therefore each reduced from the values 
given in (|60|) by a factor of p. The quantity 779 is in- 
dependent of p. Substituting these values in the Eq. H55|) 
and H56|l . one finds, in the clean limit, for 7^ 1: 



A^ 



Btt 1 -1- pN'^ 



{(1. 

p{N^ 



3P)(1-P) (63) 
1) [3p - 1 + 2p{N - C - 3)] } . 



Eq. (|63|l reproduces the "universal intrinsic value" H57|) 
for p = 0, and the small- angle scattering limit H61I) for 
p = 1. Note, however, that it does not simply interpolate 
between the two values: Eq. (|63|) is not necessarily a 
monotonic function of p. 

In general, the parameters A^ and C are related to each 
other. However, they enter Eq. (|63|) only via the combi- 
nation N — C,. In order to illustrate the behavior of our 
model for different values of N — C^, we plotted asnip) at 
different values of C, keeping A^ fixed (see Fig. Note 
that spin-Hall conductivity can even change its sign for 
N — C, > 2. Also, can exceed the "universal intrinsic 
value" 3e/87r. 

The model with a combination of isotropic and small 
angle scattering can also be considered as an approxima- 
tion to a situation where there is a single type of impu- 
rity with a scattering kernel that is neither isotropic nor 
strongly peaked at small angles. 



VII. RELATION TO EXPERIMENTS AND 
PREVIOUS THEORETICAL RESULTS. 

In agreement with previous theoretical work, we have 
found that there is no spin Hall conductivity in the case 
A^ = ±1, applicable to a 2D electron system with pure 
Rashba or k-linear Dresselhaus coupling, for arbitrary 
form of the impurity scattering, while the spin Hall con- 
ductivity has the full intrinsic value for the case A^ = 3, 
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which models a 2D hole system, when the impurity scat- 
tering is isotropic. For angle-dependent scattering, with 
iV = 3, we find that the spin Hall conductivity is altered, 
so that its precise value, and even its sign, can depend 
on such details as the ratio of small angle to large angle 
scattering, and the energy dependences of the spin-orbit 
coupling and of the hole velocity. For the 2D hole systems 
that actually occur in GaAs, an accurate computation 
should also take into account warping of the Fermi sur- 
face due to tetragonal anisotropy, which has been omitted 
from our model. 

Our results imply that spin-Hall conductivity for N ^ 
1 is not universal. Quantization of ctsH in units of e /Stt 
is broken by small-angle scattering processes even in the 
clean limit [boT 3> 1), when spin-orbit field is strong. In 
contrast, charge Hall conductivity in the strong field limit 
is insensitive to the details of disorder scattering. This 
difference probably reflects the non-topological origin of 
spin-Hall conductivity. This conclusion is indirectly sup- 
ported by the result of numerical Laughlin gauge flux 
experiment performed in Ref. |35j| for N = 1. 

In a recent article, Liu and Lei2& have performed nu- 
merical calculations of the spin Hall conductivity for a 
two-dimensional hole model, with scattering from impu- 
rities set back 50 nm from the layer. Their model has a 
quadratic energy spectrum, and k"^ spin-orbit coupling, 
corresponding to our model with N — N — 3, C, — 0, and 
has densities varying from lO^*' to IQ^^cm^^. (The den- 
sities in the published version of Ref. were misstated 
by a factor of 100). Their numerical results for CTsh range 
from slightly smaller than the intrinsic value, at their 
lowest densities, to roughly twice the intrinsic value at 
high densities, in all cases quite different from the pre- 
diction of our Eq. H62() . For high densities, it appears that 
their model is in the regime where the spin-orbit splitting 
b/vp is larger than the momentum transfer in a scatter- 
ing event, so that our formulas would not apply. (The 
numerical results actually coincide with a prediction for 
this regime in a new work by Khaetski:'^^) However, we 
do not have a complete understanding of the numerical 
results at lower densities. 

Although we do not know the precise value to be ex- 
pected for the spin Hall conductivity in the 2D hole sam- 
ples studied by Wunderlich et al.f^i it seems plausible that 
the magnitude should be similar to the universal intrinsic 
value for N = 3. In recent work, Nomura et alm^ have 
analyzed the geometry of the experiment in Ref. Q and 
conclude that the amount of polarized light obtained is 
consistent in magnitude with what might be expected to 
arise from the predicted intrinsic spin Hall conductivity. 
Nevertheless, more work seems necessary before one can 
be confident that one has a complete understanding of 
the origin of the the observed polarization in these ex- 
periments. As has been noted by various authors, spin 
polarization near a boundary depends on the boundary 
conditions, and may not be simply related to spin Hall 
currents away from the boundary. As an example, for the 
2D electron system with pure Rashba coupling, where 



fsH = 0, one may still find spin polarization in the z- 
direction near the lateral boundaries of a sample if spin- 
fiip processes are strong at the boundary)^ or if electrons 
can cross the boundary into a contact or a region with dif- 
ferent spin-orbit coupling^®-'^^. In any case, it is expected 
that the polarization near a boundary will be insensitive 
to processes that occur further from the boundary than 
a few times the length scale for spin relaxation, and this 
length scale is quite short for holes in GaAs. It seems 
unlikely that one will be able to obtain a quantitative 
measure of the bulk spin Hall conductivity in such sam- 
ples based on observations of spin accumulation near a 
boundary. 

In an earlier theoretical work by two of the present 
authors f21 a set of Boltzmann-likc kinetic equations was 
derived for a matrix distribution, denoted /p(x, t), which 
contained spin information but was integrated over the 
energy e. The model was restricted to the case of k- 
linear spin-orbit coupling, and impurities were not in- 
cluded. The focus of Ref. jH3 was on the time evolution 
of distributions that are non-uniform in space, and the 
linear response to a time-dependent perturbation in the 
absence of impurity scattering. The formalism developed 
there cannot be directly applied to spin-dependent trans- 
port in the presence of impurities, which is the focus of 
the present work. 
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APPENDIX A: DERIVATION OF SCATTERING 
KERNEL AND THE ELECTRIC FIELD TERM 

1. Scattering kernel 

We give here a derivation of the two terms M'^and M™ 
in the scattering kernel 1)31(1 , using an approach based on 
Fermi's Golden Rule. 

The term M'' arises from the spin-orbit corrections 
to the density of final states. Consider an unpolarized 
electron of energy e, traveling in a direction 9' , incident 
on an impurity located at a point Tq. (The incident 
electron may be considered to be an incoherent 50-50 
mixture of spin up and spin down states along any con- 
venient axis of quantization.) Suppose the electron is 
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scattered into a direction 9. At the angle 9, there are 
two possible momenta + 5k, with Sk = h(9)a /2v^, 
where cr = ±1 denotes spin parallel or antiparallel to 
h{9). As in the discussion following l|10|l . the corre- 
sponding densities of states, per unit 59, are given by 
(27r)^^(fcc + 5k)/{ve + 5v), and the difference between 
(7 = ±1 is given by 5v = {2Tr)~^ {d / de){k^b/Vf) . If we 
ignore possible differences in matrix elements for the two 
spin states, the difference in transition rates is given by 
W{q)5i>/2, where W{q) is defined in This leads to a 
net polarization of the scattered spin, giving a contribu- 
tion to d^{9)/dt equal to M'^{9,9')5ni9'), where M'' is 
given by H32(l . Since spin is conserved during the scatter- 
ing event, we will also have a contribution to d^{9')/dt 
equal to -M'^{9,9')5n{9'). 

The term M"" arises from spin-dependent momentum 
shifts in the matrix elements. The momentum trans- 
ferred during the scattering process depends on the spin 
orientation of both the initial and final state. Letting 
ct' = ±1 distinguish between initial spin parallel or an- 
tiparallel to h{9'), we see that the actual momentum 
transfer is g + 5q, where q is given by Eq. H25|) and 

5q = (1/4) [b{9)a + b{9')a'] (dq/de) (Al) 

is due to spin-orbit interaction. As we are assuming that 
b/Ep <^ 1, and we only need the scattering kernel correct 
to linear order in 6, we may consider the effects of b{9) and 
b{9') separately, and add the results at the end. Suppose 
that b{9) = 0. Then 5q = [b{9')a' /4:]{dq/de). The change 
in the squared matrix element for the scattering process 
is given by 5W — [dW / dq)5q, which leads to different 
scattering rates for a' = ±1. The out-scattering process 
thus creates a net polarization of the electrons remaining 
at angle 9' , and gives a contribution to d^{9') /dt) equal 
to 



''^''^■t.n('-^]'J^^5ni9') 



89 



(A2) 



Since spin is conserved in the scattering process, we must 
have a contribution to d^{9)/dt which is the negative of 
(IX2I1 . 

Next consider the case where b{9') = 0. Now 5q de- 
pends on a, so the scattering rate depends on which of 
the two final states is involved. The difference in matrix 
elements leads to a contribution to the polarization at 9 
given by 

d^{9) h{9) f9-9'\dK{9-9')^,^,^ , , 

- tan I — - — I -5n{9'). (A3) 



dt 



89 



Again, since spin is conserved in the scattering process, 
we must have a contribution to d^{9')/dt which is the 
negative of ljA3|) . The sum of these four contributions 
give the contributions proportional to M"'(0, 9')5n{9') in 
Eq. (jSH). 

Although we have used Fermi's Golden Rule and the 
Born approximation in deriving these results, we expect 
that Eqs. H32(l and (|33|l should hold more generally, 
within the model given by the Hamiltonian as long 
as the impurities are sufficiently dilute. 



2. Electric field term, Eq. 

The direct effect of a uniform electric field term in the 
Hamiltonian, in an infinitesimal time interval 5t, is to 
displace the momentum distribution by an amount (5k — 
Ee<5<, while keeping fixed the orientations of the spins. 
To see the effect of this, it is convenient to integrate over 
energy and the magnitude of the momentum, and define 
a spin density at an angle 9 by 

niW^y^ ^ j dE[a^pnp^{\,.e)l (A4) 

with k oriented in the direction 9. In thermal equilib- 
rium, m is given by 



m°^(6l) = h{9)kF/{^TT'^VF) 



(A5) 



under the assumption that 6 and the temperature T are 
small compared to Ep. 

The electric field term conserves the total spin, but 
gives rise to a flow of m around the unit circle at a ve- 
locity 9 = —Eekp^ sin9, for E along the x-axis. Then, 
to lowest order in the electric field, we have 



dni{9) 



at 



Ee d 
8n^VF 89 



{hsm9). (A6) 



By definition, [m{9) - m'"i{9)] = J de^{9,e). More- 
over, it is clear that # = for energies far from the 
Fermi energy, and for any fixed value of 9, we must have 
d^/dt\E oc —dfa/de. Equation 129|l then follows from 
"[1 



APPENDIX B: ALTERNATE DERIVATION 
USING GREEN'S FUNCTIONS 

The kinetic equation H26|l - (|33|l contains non-diagonal 
elements of the spin density matrix and is therefore a 
quantum kinetic equation, rather than a conventional 
Boltzmann equation, which contains only occupation 
numbers. In this section we show that this equation 
can be derived from the microscopic Hamiltonian Q by 
Green's function methods. In our derivation, we shall 
mainly follow the route outlined in Ref. '4l'|. As usual, 
we introduce one-particle retarded and advanced Green's 
function Gr and Ga, and Keldysh function Gk describ- 
ing non-equilibrium state of the system. (Note that all 
Green's functions are operators in spin space.) Treating 
the disorder potential as a perturbation, one may write 
Dyson equation for these function in a matrix notation: 



G^ 
G^ 



(Bl) 



Here the lower bar denotes matrix in Keldysh space, and 
Gq^ = idt — H. In this work, we neglect localization ef- 
fects, and consider only diagrams with non-crossing im- 
purity lines. This gives the standard approximation for 
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the self-energy part 

E(x, t, x', t') = W?(x - x')G(x, t, x', t') , (B2) 

where W{x — x') = (V^(x)F(x')} is the disorder correla- 
tion function. 

For a time-independent disorder and in the absence of 
electron-electron and electron phonon interactions, the 
functions Gu and Ga are independent of Keldysh func- 
tion Gk- In fact, they are Hermitian conjugate to each 
other: 



GA{E,k)^&^{E,k) , 



(B3) 



and it is sufficient to find only one of them. Since 
the functions Gr and Ga are independent of the non- 
equilibrium state, they depend only on x — x' in a trans- 
lationally invariant system. Going to the Fourier repre- 
sentation, one finds for Gr: 



Gfl(i?,k) 



1 



, (B4) 



S-ek + iT.b(k) -I]fl(i;,k) 
where ek is the hole dispersion without spin-orbit, and 



W{k-k')GB,iE,k') , (B5) 



where M^(k) is the Fourier transform of the disorder cor- 
relation function W{x — x'). 

In general, the self-energy part l|B5|l is an operator 
in spin space. Its Hermitian part describing disorder- 
induced renormalization of electron spectrum and spin- 
orbit coupling will be ignored in the following. The anti- 
Hermitian part of T,ii{E,k) which we denote —i/2f de- 
scribes the decay of one-particle state because of elas- 
tic scattering. It is convenient to project the function 
GR{k',E) onto spin eigenstates: 



GR{Ey) = 



1 



,=^E-E,{k') + ^:^. 



I3(T ■ h{k') 



where 



E^ik)=e^-^b{k) 



(B6) 



(B7) 



are corresponding eigenvalues, and r^(i?, k) are lifetimes. 
Since the imaginary part of GR{E,k') is peaked near 
E ~ Ei3(k'), one can replace it by a delta- function to get 
the self-energy in the leading order in impurity concen- 
tration. Thus, we find 



r{E,k) 



= 2tt 



(Pk' 



■El 



/3cr • b(k') 
Kk') 



(B8) 



X W{k~k!)5{E - Ep{k')) 



In this work, we restrict ourselves by the first order cor- 
rections in spin-orbit coupling. Also, we will be primarily 
interested in the decay time on the "mass shell" : 



E = e,, 



a- ■ b(k) 



(B9) 



Expanding Eq. (|B8|) to the first order in b(k'), one finds, 
on the mass shell. 



1 



27r 



(fk' 



■W{k-k!)[5{€^~e^,) (BIO) 



- *T.(b(k)-b(k'))<5'(ek-ek')] 



Thus, in general, scattering rates are different for spin 
majority and spin minority bands. 

Having found the retarded function Gr{E, k), one can 
now rewrite the remaining equation for Gk to make it 
uniform in Gk (see Ref. ^l| for details): 



G^— I/O /^^i 



GrTjk — ^kGa , 



(Bll) 



where 



SK(x,x',t,i') = M/(x-x')GK(x,x',t,t') , (B12) 

and operators G]^^ and act on the left and the right 
argument of Gk respectively. We now apply Wigner 
transform: 



G^(t,x;i',x') = / ^ 



dE d'^k 



iU.-{x-x')-iE{t-t') 



27r (27r)2 

t + t' x + x' 



X 5kB , 2 ' 2 

and the left hand side of Eq. IjBlip becomes 

{G^Gk — GrGj^) 
%E , 1 f k ld{h- a) 



(B13) 



(B14) 



dt 



dk 



\ [cr • b(k),gk£;] + M ^,5k£; 
2 2 I Tk 



where the scattering rate operator 1/fk is given by 
Eq. ||BT0)| . The right hand side of Eq. (|BTT|) takes the 
form 



(GrX'K — ^kGa) 
d^k! 



(B15) 



W{k-k') GR{E,k)g^,EM 



gw'Ei^,t)GA{E, k) 



(here we assumed that x- and t- dependence of gkB is 
smooth on the scale of Fermi wavelength) . 

Calculating transport quantities, we are normally in- 
terested in Keldysh Green's function at coinciding tem- 
poral points, t = t' , which in Wigner representation 
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corresponds to energy-integrated functions. However, 
Eqs. HBll() - (|ijl5|) contains Green's function for arbitrary 
energies. To express all quantities in terms of functions 
at t = t' , we use the following ansatz for g]iE- 



gkE^GRiE,k)hk~h^,GA{E,k) , 



(B16) 



where the function depends only on momentum. In- 
tegrating Eq. (|B16|I over E, one finds: 



dE ^ J 
Ztt 



(B17) 



Therefore, the function hk{x,t) can be used to compute 
transport quantities. In fact, the function /k = 1 — 2/ik 
can be interpreted as a distribution function, with diag- 
onal elements of /k being the occupation numbers. 

Now we can derive the equation for the function 
/ik(x, t) by integrating left- and right-hand sides of quan- 
tum kinetic equation ljB14() and ljB15|l . Since the left- 
hand side depend on E only via Gk{E), the integration 
is trivial. The terms on the right hand side (Eq. (jB15|) ) 
contain E explicitly, and should be carefully integrated. 
Substituting the definition l|B16|l into the Eq. llB15|l . one 
gets four different terms. Two of them contain the prod- 
uct of two retarded (or two advanced) Green's functions, 
and vanish after integration over E. The two remaining 
terms are anti-Hermitian conjugate to each other. One 
of them has the form 



Ri 



—GHiE,k)hk,GAiE,k') 
Zn 



(B18) 



Projecting Green's function onto eigenstates with the 
help of Eq. ljB6|) and integrating over E, one finds: 



/3,/3'=± 



^cr • b(k) 



(B19) 



!+/?'■ 



Ep,{k') 

T-h(k') 
~b(k^ 



r/j(k) 



'(k' 



where the spin-majority and spin-minority energy bands 
are given by the Eq. ljB7|l . Assuming the scattering rate 
I/tct in the denominator to be small, we now rewrite this 
term using Sokhotsky formula: 



1 1 ■ N 

-— = P - - ZTTd X 

+ lO X 



(B20) 



where P denotes principal value. Subtracting Hermitian 
conjugate and expanding to the first order in spin-orbit 
field b(k), one finds, for the right hand side of quantum 
kinetic equation: 



Ri-R[ 



27r/ik'(5(ek - fk') 
Tr[hk',(T-{h{k)^ 

\',T-(b(k) 



iP 



h{k'))}s'{ 
fb(k')) 



Ck 



(B21) 



(ek - Ck' 



The first term in the Eq. (|B2ip is a conventional scatter- 
ing term, while the second is the spin-orbit correction to 
the scattering rate. The third term is due to renormaliza- 
tion of one-particle spectrum by the disorder, and should 
be ignored in the following, since we already neglected 
similar terms in the real part of self-energy. Now we can 
use the Eq. ljB21|l and ljB14|l to write the equation for 



dt 

i 



1 f k la(b-cr) 
dk 



2 

2n 



2 [m 2 
cr • b(k), ft-k 
d^k' 



(B22) 



W{k-k')[{hu' ~hu)S{e'^ 



Ck' 



- \ {/ik' - /ik, T • (b(k) - b(k'))} 5'(ek - ekOj 

(we grouped together terms proportional to W). 

To this point, we did not take into account the effects 
of external electric field. To incorporate the electric field 
into kinetic equation, consider a very smooth external 
scalar potential U{r). Strictly speaking, this potential 
breaks translational invariance, and the calculations we 
have done so far become invalid. However, for small field 
the potential U(r) is smooth on the spin-orbit precession 
length hve/b and does not lead to transitions between 
spin subbands. Therefore, U{r) can be considered as a 
local shift of electron energy. The term G^^ G k — G kG^^ 
now contains U{r) — C/(r'), which becomes, after Wigner 
transform, VrC/VkgkB- Thus, the electric field term in 
the kinetic equation takes the standard form. Rewriting 
kinetic equations in terms of function /k = 1 — 2/ik, one 
finds: 



dt 



k ld{h- cr) 



9k 



V/k 



(B23) 



T-b(k),/k 
d^k' 



eE 



dk 



I {/k' - /k,<T • (b(k) - b(k'))} (5'(ek - ek') 



The kinetic equation ljB23p is sufficient to compute spin 
transport in a 2D hole gas. Note, however, that in gen- 
eral, spin-orbit terms on the right hand side mix states 
with different values of Ck because of the presence of the 
derivative of delta-function. Thus, in general, one has to 
consider all electron states with different energies simul- 
taneously. This mixing occurs because the quantity ek 
is not conserved in the scattering process. Instead, the 
true electron energy 



Ek 



b(k) 



(B24) 



is conserved. We are now going to use this integral of mo- 
tion to separate contribution of electrons with different 
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energies, and make energy conservation explicit. We shall 
make the change of variables, expressing the electron mo- 
mentum in terms of energy e and momentum direction 
9. Note, however, that the electron energy e contains 
spin-dependent part, and therefore replacing momentum 
length by e is not a simple algebraic transformation: one 
has to ensure the proper ordering of spin operators in the 
resulting expression. To figure out this operator order, 
consider the effect of spin-dependent gauge transforma- 
tion 

^'{r) = Ua'iiir) = exp(i(q • r)(cr • b))*(r) (B25) 

on the one-particle density matrix p^.r' = 5'^(r)\E'(r'). In 
the first order in b, one finds: 



iq • (r — r') 



iq- (r + r') 

Pr,r ^ [Pr,r' , CT • bj 



{/5r,r', cr • b} 



(B26) 



To find the electric field contribution, note that for small 
field it should be computed only for an equilibrium dis- 
tribution. In the equilibrium, the occupation numbers /k 
depend only on the total energy. Thus, 



dl 
dt 



eE^^^eE^[v,{K)cose (B32) 

a ■ h{k^) dv^ V, d{h ■ cr) 

■ cos — cos 



2 de 
I9(b -cr) sin 61 

2 



de 



do 



(the second term here takes into account the momentum 
shift given by Eq. (|B29|) ). 

In the scattering term, the delta-function and its 
derivative combine into (5(e — e'), due to energy conser- 
vation. In zeroth order, the scattering contribution is 



By doing Wigner transform with respect to r and r', one 
arrives at: 

p'(R,p) - /5(R,p)-i[/5(R,p),(q-R)(^T-b)] 

+ {(q-Vp)p(R,p),(*T-b)} . (B27) 



The second term in the Eq. (|B27|) describes spin rota- 
tion due to transformation Ua, and the third term de- 
scribes spin-dependent momentum shift. Note that spin 
operator in the last term acts symmetrically on p. This 
symmetry is related to the antisymmetric action of mo- 
mentum operator on ^(r) and ^'^(r') in Pr,r'- One can 
undo the spin rotation described by the second term in 
the Eq. ljB27p . and use the third term as a prescription 
of how the spin-dependent momentum shift should be 
performed: 

/5(R, p + q(crb)) = /3(R, p) + i {qVp/5(R, p), (crb)} . 

(B28) 
I by 

(B29) 



Now we can replace momentum in the Eq. ljB23ll by 

crb 



kfk - 



2v, 



■k , 



where k = (cos 0, sin 0) is the unit vector in the momen- 
tum direction, and k^^ and v^^ — de^/dk are defined in 
Sec. mil To rewrite the gradient term in this represen- 
tation, one can simply replace the momentum according 
to the Eq. (|B29|) : 



dl 
dt 



2V'''^ 2 ^ 2 de ' ^"-^ 



(B30) 

The spin-orbit rotation term already contains spin-orbit 
field and should not be modified at all: 



9/ 
dt 



T-b(k),/ 



(B31) 



dl 
dt 



k,de' 

2TrVe 



ly(k-k') (/>')-/>) 



(B33) 



where both k and k' are taken on the surface — e. 
Spin-orbit corrections to the Eq. ljB33|l arise from either 
volume element (Pk.' or the matrix element W^k — k'). 
The matrix element contribution is due to spin-orbit in- 
duced changes of k and k' : 



dl 
dt 



k,de' 



(B34) 



Transforming this expression as explained in Appendix 
A, one recovers the scattering kernel M^u. 
The volume element is computed as follows: 



d'k = kdkde 



K + ^) fd9 . (B35) 
zVf I de 



The derivative dk/de should be computed to the first 
order in spin-orbit field, as explained in comments to 
Eq. UlUf) . and the volume element becomes 



Vf 2 de Vf 



(B36) 



(as usual, its matrix part should be applied symmetri- 
cally). This gives additional contribution to the collision 
integral: 



dl 
dt 



d k,(T ■ b(k^) 



(B37) 



de 2ve 



J{e,9')-f{e, 
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The evolution of the density matrix p in e ~ represen- 
tation is then given by 



dh 


df 


dt 

a 






Of 




df 


'dt ~ 


dt 


at 

b 


^ 

E 


dt 


^ 




dt 



dl 
dt 



(B38) 
fR32ll . 



where ah terms are given by Eqs. I|B30() . (|B31() 
(fR34ll and (ImTll . 

To simpHfy calculation of physical quantities, one can 
remove the spin dependence from the volume element. 
To do that, we include the volume element in the distri- 
bution function via the following transformation: 



n(e,0) = ^/(e,( 



d k^cr ■ h(kf) 

de Vf 



(B39) 



One can see that in the first order in spin-orbit 

^ {fie, e),d^kj = n(e, e)dede . (B40) 

The function n(e, 9) can be now identified with the 
distribution used throughout the paper. Substituting 
Eqs. (|B39l) and ^ into Eqs. p32l) . (|B3l and p37|l . 
one can reproduce Eqs. 129() . (|33l) . and H32() respectively. 
Indeed, the spin-orbit correction to the electric field term 
is given by the sum of contributions from Eq. ljB32p and 
Eq. iB39t : 



dn 
'dt 



-eE 



dfa (T 

de 2 



Vf cos 9 



d kfhf 



kfhf dve 
cos9—- 

Vf de 

sin 9 dh 



kf cos 9 



de Vf 
dhf 



(B41) 



de 



d9 



Collecting the terms proportional to cos9, one finds: 



dn 
'dt 



E,so 



r^dfo a- 
^^572 



cos 9 — 

Vf 



^9dh 



d9 



(B42) 



thus reproducing Eq. H29|l . The matrix element correc- 
tion to the coUision integral (|B34|) reproduces Eq. 
The density of states correction (|B37|l . however, is modi- 
fied by a similar contribution that arises from transforma- 
tion of the collision integral ljB33|l . In the latter, one has 
to express /(e, 9') in terms of n{e, 9) and then transform 
the resulting contribution according to Eq. ljB39|l to get 
dh/dt. As a result, four spin-orbit terms arise. Then, 
the term in Eq. (|B37I) containing fi{9') is cancelled by 
the corresponding contribution to the Eq. l|B33p . Similar 
cancellation occur for n{9) term. As a result, one re- 
covers Eq. H32|) . This concludes the derivation of kinetic 
equation H26|l by Green's function method. 
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